Phase heterogeneities of lipidic aggregates 
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We propose a model for explanation the "domain- wall" type configuration states in binary lipid 
mixtures of cationic and neutral lipids, associated with observed relaxation effects in their aggre- 
gates. We apply the analogy with formation of Kibble-Zurek topological defects, which we suppose 
connected with structural dynamics of the lipid phases. In frames of the proposed model, the density 
of kink-type defects and the energy of the configurations are calculated. 

PACS numbers: 



I. INTRODUCTION 

A number of applications for lipids and their mixtures 
in the biomaterial technology, therapy and industry cause 
the necessity of careful theoretical predictions of their 
phase transformations. Lipidic aggregates, as well as ly- 
otropic systems at all, are ideal substances demonstrating 
the quantum topological phase transitions ( [H, H| and ref- 
erences therein). Some experiments on polar fluids 0-0| 
reveal new interesting phenomena of self-organization. 

For concentrated suspensions, pastes, emulsions, 
foams, and associative polymers, the mechanism of struc- 
ture relaxation in soft solids, which is based on the me- 
chanical anelastic spectroscopy in rheological frequencies, 
was proposed by Wyss and coworkers [3(. Amplitude- 
dependent measurements have shown, that when the 
strain rate becomes large, it can itself drive the slow 
structural relaxation process at the time scale of the im- 
posed strain rate Also, under an applied strain, the 
observations of high-frequency shifts cause to interpret 
the enigmatical slow relaxation dynamics in a principally 
new way. While, the authors [3| seek an explanation of 
their measurements in analogous phenomena in super- 
cooled fluids 

From the other side, the recent anelastic spec- 
troscopy studies Q of dynamical processes on neutral 
DMPC (dimyristoil phosphatidylcholine), DOPE (di- 
oleoyl phosphotiylethanolamine) lipids, cationic DDAB 
(dimethyldioctadecylammonium) , DOTAP (di-oleoyl 
trimethylammonium) lipids, and their DOTAP/DOPE 
and DOTAP/DMPC mixtures brought out the hypothe- 
sis of new micro- and nanoscale structure heterogeneities 
in the lipid membranes which look like onto domain 
walls. These observations have been carried out in a wide 
temperature range starting with the supercooled state, 
and at wide-range excitation frequencies (10 2 -10 4 Hz), 
with the lipids deposited on a solid substrate. 

Weak frequency-dependent shifts in the elastic modu- 
lus are observed at these cryogenic temperatures; peaks 
on the relaxation curves at the low-dynamics regime are 
closed by the typical "smooth" relaxation, so, the ob- 
servations evidence collective short-range motions of the 
lipids [5|. 

Also, one may indirectly compare these data with neu- 



tron scattering ((J and atomic-force microscopy (AFM) 
on some neutral and cationic lipids and their mix- 
tures. 

Hence, a number of soft matter phenomena can be de- 
scribed in the common interdisciplinary modeling associ- 
ated with a mechanism of forming of topological defects, 
such as domain walls, or vortex strings. 



II. FORMATION OF TOPOLOGICAL DEFECTS 

Generally, in the cases specified above, the phase tran- 
sitions are continuous. The Kibble-Zurek (KZ) mecha- 
nism of formation of the domain-wall type defects in an 
adiabatic regime is convenient for a description of dy- 
namics of these transitions 0, Q . 

As topological defects, domain walls are broadly known 
in their universality [ToL Hlj ]. from cosmological models 
to the theory of condensed matter [l], Ojj]. The initial 
Landau-Zener (LZ) Hamiltonian [l3[ of a two-level sys- 
tem expresses the dynamical processes at classical phase 
transitions. Afterward, the LZ Hamiltonian has been 
generalized onto quantum phase transitions. 

Here, we identify the domain-wall defect formation 
with the some quantum phase transition at the last (" adi- 
abatic") stage of the LZ-type evolution, which includes 
three regimes in general Q. Then, one may apply a for- 
malism of adiabatic quantum computations, or Quantum 
Annealing (AQC-QA), with the KZ approach [12[ in ab- 
sence of frustrations. 

This modeling enables us to estimate a density of 
kinks and a residual energy, corresponding to the 
one-dimensional quantum Ising system with the time- 
dependent term of a transverse field, the Hamiltonian 
of which is [12| : 
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here, erf' 2 are Pauli matrices for the i-th spin of the chain. 
Ji denote random couplings between neighboring spins, 
and hi is a random transverse field. The function of 
time T(t) serves for rescaling a transverse field hi at an 
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annealing rate r : 

r(*) 
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One should note, that following the KZ scenario, in the 
end of the ordering into a non-equilibrium state, the tran- 
sition time tq and the average finite ordered domain size 
are connected by 



(3) 



where z and v are critical exponents [T3 |. 

The LZ Hamiltonian describes the evolution of a sys- 
tem in time t 
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In the adiabatic-impulse approach, with the evolution 
of a system in time t from —00 to 0, after completion of 
the transition, the density of kinks equals: 



n = lim 

N=+a 
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(5) 
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In other words, dynamics in a system can be exactly 
described by a series of uncoupled LZ systems [HI, . 

The authors of [HJ have analyzed the fermionic Hamil- 
tonian received by applying the Jordan- Wigner transfor- 
mation to the Hamiltonian (1). To treat the question of 
dynamics, they have solved a system of linear differen- 
tial equations received from the fermionic Bogoliubov's 
equations b y m eans of the known ansatz. So, in agree- 
ment with [14], at a finite r and at LZ factor, scaled 
by Bogoliubov-de Gennes transformation, the density of 
kinks and the residual energy, reduced to the total lat- 
tice size L, can be estimated respectively by the next 
relations [12j: 
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The C parameter has been found numerically ~ 3.4 ± 
0.2 [12J. 7 are the Bogoliubov operators diagonalizing 
H(t), n is the universal function, and II -1 denotes the 
inverse function of II [l2[ ; the characteristics g is defined 
by the equality g = — log(— Ai)/-\AL, Ai = 2 (ex + £2) is 
the excitation energy of single-particle eigenvalues t\ < 
e 2 < ■■■eh- E res = E t - E c iassi caU here E t denotes a time- 
evolved state energy, and L e (r) is a length of the defect- 
free region upon annealing. The critical point probability 
(ibid) is 
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FIG. 1: Average residual energy per site as a function of the 
annealing rate r. 



g c denotes the characteristics g in the critical point. The 
exact chain's P depending on concentration is known 

M) 

In principle, we know the classical 2D Ising simulation 
with the Glauber dynamics @ (the heat bath algorithm) 
for a non-equilibrium system under the KZ mechanism. 
The continual version of the Hamiltonian with pure re- 
laxation time is given there. It seems useful for our goal, 
because, in such model, the "domain walls" are always 
annihilating Q. However, we can not follow it directly 
by virtue of the reasons shown below. 



III. NUMERICAL MODELING AND RESULTS 



According to the hypothesis of displacing lipids motion 
Q, we carry out the numerical experiments in the spirit 
of the quantum model of [l2j and references therein, but 
for a 3D Ising lattice allowing frustrations. 

To involve "concentration" in this modeling, we have 
to keep a number of particles during simulations. 

At free field parameters, let us assume that the Hamil- 
tonian (1) is a bosonized Hamiltonian of our particles, so 
that we operate with the spin variables cry = ±1. 

Then, above the critical point [T^ |. at the periodical 
boundary conditions on the 6cc-lattice of 48 x 48 x 48 
and 60 x 60 x 60 sizes, we calculate this energy (Fig.l) 
and the kink density (Fig. 2.). 

For the Hamiltonian (1), we have performed a typical 
Monte Carlo (MC) algorithm, in which the probability 
of states is exp(— AH/ksT), where T is absolute tem- 
perature, ks is the Boltzmann factor, and AH is the 
energy difference between neighboring states. A number 
of MC running steps equals 10 6 , and the thermalization 
is carried out in 10 5 steps. 
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FIG. 2: Density of kinks as a function of the annealing rate 
r for the 60 x 60 x 60 lattice. 



IV. DISCUSSION 

At the adiabatic regime (Fig. 1), the average resid- 
ual energy behaves in a gree ment with the KZ mecha- 
nism (see, for instance, [12|]). These data allow us to 
discuss, how far the LZ theory is satisfied to the hypoth- 
esis, following from the experiment where a density 
of kinks could be measured hereafter. The lattice mod- 
els for lipids are widely known (l8j ; in their frames, an 
aggregation scale and a lattice size are comparable. 

Thus, it will be possible to construct a structure pa- 
rameter of fluid aggregates and lipidic mixtures in terms 
of quantum Ising models Q . 

Some interesting 3D calculations [19| are known for 
vortex strings in the He isotopes, which qualitatively 



agree with our results (Fig. 2.). 

However, our approach encloses the next problem. The 
KZ-type models contain a concentration dependence, an 
absence of which was emphasized especially for the cryo- 
genic experiments Q on the lipid mixtures. Sometimes 

pi El, 

at the calculations of domain wall sizes in the 
one-dimensional case, this question is imperceptible, but 
is not solved for different types of soft solids [3j, [5| . 



V. CONCLUSION 

So, in contrast to lipidic phases without defects, where 
phase transformations may be characterized in classical 
Ising models [2(| , the case of unexplored low-temperature 
phase transitions compels us to involve the quantum lat- 
tice model with a random transverse field. 

If the domain walls are annihilating and/or generating 
a new phase, then it is reasonable to continue the exper- 
iments in a wide amplitude and frequency range. The 
discussed model can be tested also in neutron diffraction 
experiments on similar lipidic membranes @. 

In the case of AFM, and in general in the presence of 
a substrate 0, [H| , the observable relaxation peaks have 
to be separately specified. 
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